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ABSTRACT 

Clusters of galaxies are enriched with positrons from jets of active galactic 
nuclei (AGNs) or from the interaction of cosmic-rays with the intracluster gas. 
We follow the cooling of these positrons and show that their eventual annihilation 
with cluster electrons yields a narrow annihilation line. Unlike annihilation in the 
interstellar medium of galaxies, the line produced in clusters is not smeared by 
three-photon decay of positronium, because positronium formation is suppressed 
at the high (> 1 keV) temperature of the cluster electrons. We show that if AGN 
jets are composed of e"'"e~ pairs, then the annihilation line from rich clusters 
within a distance of 100 Mpc might be detectable with future space missions, 
such as INTEGRAL or EXIST. 

Subject headings: galaxies: jets - X-rays: galaxies: clusters - galaxies: clusters 
- cosmic rays 

1. Introduction 

Despite decades of intense study, the composition of relativistic radio jets remains enig- 
matic. While the existence of synchrotron-emitting electrons is secure in both quasars (Begel- 
man, Blandford, & Rees 1984) and microquasars (Mirabel 2001), there is no conclusive 
evidence that can determine whether positrons or protons make up the positively charged 
component. Naively, one might expect these cases to be easily distinguishable through an 
observational search for the electron-positron annihilation line. However, modeling of the 
annihilation process in active jets has shown this not to be the case: the annihilation spec- 
tral feature is not a line but is instead very broad and hence difficult to unambiguously 
identify (although it may possibly contribute to the observed 7-ray spectrum of blazars; see 
Boettcher & Schlickeiser 1996). Even if positron cooling to non-relativistic energies within 
the jet is efficient, the high bulk Lorentz factor (7jet ~ 10) of the jet would likely smear 
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out an annihilation line (Begelman, Blandford, & Rees 1984). The resulting radiation may 
be detectable with upcoming space observatories, although its interpretation will depend 
critically on proper modeling of the region inside the jet (Wang, Durouchoux, & Li 2001). 

The chief obstacle to producing a recognizable annihilation-line feature is the highly 
relativistic nature of the associated plasma. This obstacle may be overcome simply by 
waiting for the active galactic nucleus (AGN) to become dormant. Once the central engine 
disappears, the material in the jet will presumably mix with the ambient medium and cool 
to the ambient temperature. In this paper, we calculate the annihilation signal from "relic" 
positrons produced by an AGN embedded in a galaxy cluster, as the positrons thcrmalize 
with the electrons in the intracluster medium (ICM). We will show that under reasonable 
assumptions an annihilation line from nearby clusters would be observable in the near future. 

A similar suggestion of thermal annihilation after escape from the accelerating source 
was made by Maciolek-Niedzwiecki, Zdziarski, & Coppi (1995) in the context of small-scale 
jets from stellar-mass black holes. However, in galactic environments with typical temper- 
atures T < 10^ K, formation of an annihilation hue at 511 keV is inhibited by the rapid 
formation of positronium, whose primary annihilation channel yields three photons. For 
example, the observations of Kinzer et al. (2001) show that 93 ± 4% of pair annihilations 
in the Galactic center occur through the positronium channel and that > 70% of the total 
annihilation energy is emitted in a broad continuum rather than in the line. We will argue 
that galaxy clusters arc ideally suited to producing annihilation lines because the character- 
istic temperature of the ICM is larger than the binding energy of positronium; hence, nearly 
all annihilations produce photons near 511 keV. 

Another important open question involves the dynamical significance of cosmic rays 
to the ICM. Faraday rotation studies have revealed that cluster cores have magnetic fields 
> 1 /iC (Kim et al. 1990; Clarke, Kronberg, & Bohringer 2001), while observations of excess 
emission in the ultraviolet and hard X-ray bands from the ICM have been used to infer 
the existence of widespread cluster fields at slightly lower levels ~ 0.1-1 /xG (e.g., Rephaeh, 
Gruber, & Blanco 1999). The pervasiveness of magnetic fields in clusters indicates that 
coUisionless shocks generated by accretion or merger events may be efficient particle acceler- 
ators (Blandford & Eichler 1987; Colafrancesco & Blasi 1998). Observations of synchrotron 
radio halos indicate that the acceleration of electrons by cluster merger shocks is a common 
occurrence (Kempner & Sarazin 2001). However, these observations probe only the electron 
component. In order to understand the dynamical effects of the cosmic rays on the cluster, 
we are most interested in the accelerated protons, both because shock acceleration may in- 
ject more energy into this component than into the electron component (Fields et al. 2001; 
Butt et al. 2001) and because cosmic ray protons do not rapidly lose their energy through 
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radiative cooling. 

Recent cosmological simulations suggest that the cosmic ray pressure might be as large 
as ~ 10-40% of the thermal gas pressure in clusters (Miniati et al. 2001b). To date, two 
diagnostics of the cosmic ray proton content have been proposed, both relying on the decay 
of pions produced in collisions between cosmic ray and thermal protons, namely synchrotron 

and inverse Compton emission from secondary electrons and positrons produced in charged 
pion decay (Blasi & Colafrancesco 1999; Dolag & EnBlin 2000; Miniati et al. 2001a) and 7- 
rays produced in neutral pion decay (Colafrancesco & Blasi 1998; Miniati et al. 2001b). The 
former process suffers from the possibility of contamination by newly accelerated electrons, 
while the 7-ray signal will remain below the detection limits at least until the launch of the 
GLAST mission^ in 2005, largely because the energy from the decaying pions is distributed 
over a very wide range of photon energies. 

Because of the rapid cooling of positrons produced in the decay of tt"*" to non-relativistic 
temperatures, we would expect an annihilation line to be produced in this case as well. (Here 
too the formation of positronium is inhibited by the relatively large cluster temperatures.) 
We therefore also calculate the annihilation spectrum produced through secondary positron 
production by cosmic ray protons. However, we find that this signal is well below the 
detection thresholds of upcoming instruments. Hence, AGN are the only realistic pollutant 
of substantial amounts of positrons into galaxy clusters. We thus argue that the future 
detection of positron annihilation line radiation from clusters would constitute a robust 
signature of electron-positron jets. 

We begin by describing the factors determining the evolution of the positron popula- 
tion in §2, including cooling, annihilation, and source terms. We then solve the evolution 
equations in §3. In §4 we calculate the resulting positron annihilation signals. Finally, we 
conclude in §5 with a discussion of our results and prospects for future observations. 

Throughout the paper, we assume a ACDM cosmology with Qq — 0.3, fl^ — 0.7, and 
Ho = 70 km s"^ Mpc'^ 

2. Evolution of the Positron Distribution P\inction 

The time evolution of A^+(7, t), the mean differential number density of positrons with 
Lorentz factor between (7, 7 -|- ^7) in the cluster at time t, is described by the analog of the 
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Boltzmann equation with source and sink terms (Sarazin 1999): 



dN+{j,t) _ d 
dt 



[6(7, t) 7V+(7, t)] + g(7) - A^+(7, t)A{^)- 



(1) 



Here, 6(7, t) = (d^/dt) is the coohng rate of a single positron, Q{'~f)d'~f is the production rate 
of positrons per unit volume in the interval (7,7 + ^7), and ^(7)^7 is the annihilation rate 
of positrons in the same interval. We have neglected diffusion and loss of positrons between 
volume elements in the ICM (Ginzburg & Syrovatskii 1964). This should be a good approxi- 
mation provided that the sources of positrons are distributed over a sufficiently large volume, 
as would be expected for injection by both AGNs and cluster shocks. Alternatively, because 
clusters trap all but the very highest energy cosmic rays over a Hubble time (Berezinsky, 
Blasi, & Ptuskin 1997), our evolution equation can be viewed as simply describing the total 
positron population in the cluster, with the appropriate averages of the cooling, source, and 
annihilation terms (Sarazin 1999). 

We arc concerned primarily with those positrons that cool sufficiently to contribute 
to line emission upon annihilation. Therefore, we wish to track the number of positrons 
that have thermalized with the ambient electrons. In the test particle approximation, a 
positron cools until its kinetic energy ex ^ O.dSksTe, where ks is Boltzmann's constant and 
Te is the temperature of the ambient thermal electrons (Trubnikov 1965). In actuality, the 
finite population of cooled positrons thermalizes with the ambient electrons on a timescale 
(Trubnikov 1965) 



where Ue is the electron number density in the cluster core and Tkev = (/cfiTe/keV). Because 
Ttherm IS much shortcr than the other relevant timescales, we assume that thermalization 
occurs instantaneously. The pool of cooled positrons essentially acts as an absorbing wall in 
the cooling equation, so we write the total positron distribution function Nt{j) as 



where 7eq is the mean Lorentz factor of the ambient electrons and n+ is the number density 
of the thermalized positrons. 

In the following subsections, we describe each of the terms on the right-hand side of 
equation (1) in detail. We solve for the N^{j,t) appropriate to our two models in §3. 





NT{^,t) = n+{t)5{^ - 7eq) + A^+(7,i), 



(3) 
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2.1. The Loss Rate 6(7, 



The principal cooling mechanisms for cosmic ray electrons or positrons in a cluster are 
synchrotron, inverse Compton (IC), and Coulomb cooling (Rephaeli 1979). The synchrotron 
and IC loss rates for a single positron with Lorentz factor are given by (Rephaeli 

1979) 



c?7 
'dt 



syn,IC 



{ub + ucmbW, 



(4) 



where (Tt is the Thomson cross section, mg is the electron mass, ub — B^/Stt is the energy 
density in the cluster magnetic field, and ucmb is the energy density of the cosmic microwave 
background (CMB). The loss rate due to these two processes is therefore 
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where z is the cosmological redshift at cosmic time t. 

The loss rate from Coulomb cooling is approximately (Sarazin & Kempner 2000) 



d7\ 



inn^e'^ ^ f 1.12m ec'^-f^^^P^ 



In 



(6) 



where (3 — {1 — 7"^)^'^^ is the positron speed normahzed to the speed of hght and Up is the 
plasma frequency of the ambient electron gas. The corresponding choice for the Coulomb 
logarithm is valid for a -C /3 -C 1, where a is the fine structure constant, and for 7 1, 
with errors < 10% in the intermediate regime (Sarazin & Kempner 2000). Evaluating the 
constants, we find 
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(7) 



where Uc 



(ne/cm 



The total loss rate is then 6(7, t) = &coui(7) + ^syn.ic(7- 1). Comparing equations (5) and 
(7), we see that IC/synchrotron cooling dominates for 7 > 200(ne/10~'^ cm~^)^/^; hence our 
assumption that 7 3> 1 in equation (4) is appropriate. 
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2.2. The Annihilation Rate ^(7) 

The annihilation rate A (7+) of positrons at a given Lorentz factor through the process 
e+e" — >• 27 may be written as (Svensson 1982) 

yl(7+) = r d^-N.{^.) f ^jln^2(3Me{lcm), (8) 
Ji J-1 ^ 7+7- 

assuming that the distribution function of electrons, A^_(7_), is isotropic but otherwise 
arbitrary. Here 7+ (7_) is the Lorentz factor of the positron (electron) in the cluster rest- 
frame, 7cm is the Lorentz factor of the positron in the center-of-mass (CM) frame, /i is 
the interaction angle in the cluster rest-frame, and o"ee(7cm.) is the total annihilation cross 
section in the CM frame. We find that assuming that the ambient electrons are all at rest 
[A^_(7_) = ne5(7_ — 1)] makes only a small difference to the total annihilation cross-section 
for the non-thermalized positrons, because realistic cluster temperatures arc highly non- 
relativistic (so that 7+ 3> 7_). In this "cosmic ray approximation" 7^^ = (7+ -|- l)/2 and 
the annihilation rate becomes 

where is the classical electron radius. This expression is valid provided Coulomb correc- 
tions may be neglected, i.e. as long as /?+ ^ a; this condition is satisfied for all 7+ > 7eq at 
the typical cluster temperatures of ksTc > 1 kcV. 

In order to calculate the spectrum of annihilation photons, we are also interested in the 
derivative of the annihilation rate with respect to photon energy. This may be written in 
the cosmic ray approximation as (Svensson 1982) 

^^=n.fd,,NMiv%(k.,,)), (10) 
where k = e^/meC^ is the normalized photon energy and 



iJ'^lk-. n '""'i f -(3 + 7+)/(l + 7+) + (3 + l+)/k - lie , 



is the angle-averaged emissivity per electron-positron pair. 

For the thermalized positrons, the cosmic ray approximation is invalid because 7+ ~ 7_. 
We therefore must evaluate the general expression for the photon production rate, 

dn^ , . r°° , „ . da. 



dk dt 



/OO POO 7 

dj+N{j+) dj.N_{j_){v—{k,j+,l-)), (12) 
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for Maxwell-Boltzmann electron and positron distributions at a temperature T^. Svensson, 
Larsson, & Poutanen (1996) show that the result is 



dk dt 



line 
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(13) 



in the approximation /c -C 1 and ksTg <S nieC?- Here we use the subscript "line" because 
annihilation of the thermalized component yields a well-defined spectral line (see §4 below). 
The total annihilation rate of the thermalized positrons is then simply 



„ dn^ 
dk 



dk dt 



8 X 10~^V s-\ 



(14) 



line 



with only a very weak dependence on temperature in the range of interest. (The factor of 
two in the denominator appears because each annihilating positron yields two photons.) 

Note that in many astrophysical situations, positron annihilation occurs principally 
through the formation and subsequent decay of positronium. In such situations, the domi- 
nant annihilation mechanism is the decay of orthopositronium, which must emit at least three 
photons in order to conserve angular momentum (Berestetskii, Lifshitz, & Pitaevskii 1982). 
This mechanism therefore precludes the formation of an annihilation line. Fortunately, in 
our case positronium formation is insignificant. For temperatures ksT^ ^ 1 keV, the cross 
section for positronium formation amounts to < 10% of the cross section for direct anni- 
hilation, and the positronium formation rate decreases strongly at yet higher temperatures 
(Crannell, Joyce, Ramaty, & Wcrntz 1976). We therefore neglect this annihilation channel 
in our calculations. Note, however, that in cool galactic environments (T < 10^ K), annihi- 
lation through positronium formation dominates and the annihilation line is suppressed. 



2.3. The Source Function (5(7) 

In the case of positron injection purely via an AGN outfiow, the source function Q(7) 
is nonzero only at the instant of injection. In this section we describe positron production 
mechanisms in galaxy clusters over and above direct injection from AGN. The dominant 
mechanism, the proton-proton channel, is described in §2.3.1. For completeness, we examine 
other mechanisms in §2.3.2. 
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2.3.1. Positron Production from Proton- Proton Collisions 

Collisions between thermal and cosmic ray protons produce positrons primarily through 
the reaction chain 

p + p n+ + X, 
7r+ + i/^, 

provided that the cosmic ray proton has a Lorentz factor exceeding 7th = 1-3, the threshold 
for pion production (Mannheim & Schlickeiser 1994). The resulting positron production rate 
is then (Moskalenko & Strong 1998) 

POO rt^'^^ 
Qpp{l)^nH d-fpJp{-fp)a^+{-ip) j d'y^F+{-f,-f^)F^{-f^,jp), (15) 

ip ITT 

where uh — is the proton density in the cluster core, Jpi'jp) is the differential cosmic 
ray proton flux at Lorentz factor jp, a^^+^jp) is the cross-section for tt"*" production in the 
collision (including all channels), F^{'y,^^) is the normalized positron distribution function 
at positron Lorentz factor 7 for positrons produced by the decay of a pion with Lorentz 
factor 7,^, and ^7^(7?^ 7p) is the normalized pion distribution function produced by a collision 
between a thermal proton and a cosmic ray proton with Lorentz factor 7p. The maximum 
cosmic ray energy 7™^^ is determined below. The other integration limits are determined by 
kinematics: 7™"^ is the minimum pion Lorentz factor needed to produce a positron with 7, 
^max -g ^jjg maximum pion Lorentz factor that can be produced by a proton with 7^, and 
7™™ is the minimum proton Lorentz factor required to produce a positron with 7. 

Our calculation of the resulting distribution function follows closely the method out- 
hned in Appendices B and C of Moskalenko & Strong (1998), so we only summarize it 
here. Unfortunately, there is no flrm theoretical understanding of pion production in proton- 
proton collisions. The production cross-sections are experimentally well-determined (see the 
compilation of Dermer 1986), but the resulting F^(7^,7p) are not. Instead, we must use 
approximate models of the interaction. At low energies (7^ < 3), the "isobaric model" of 
Stecker (1970), in which the collision forms a A isobar which subsequently decays into a 
pion, fits the data best, while at high energies (7p > 7), the "scaling model" of Badhwar, 
Golden, & Stephens (1977) provides a better fit. In the intermediate regime, we take a linear 
interpolation between the two models. These models do not include the deuterium channel 
ipP ~^ D7r~^), important for low-energy collisions, and we include this separately (note that 
this channel can be treated exactly). We assume that are produced isotropically in the 
TT'^ decay, and we include polarization in the subsequent /i'^ decay (see Moskalenko & Strong 
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1998 for details). Note that F+(7, 7^) can be found exactly, unlike the pion distribution 
function. 

Our treatment neglects two effects in pion production. First, we ignore the kaon channel, 
in which the collision forms a that subsequently decays into either a tt"^ or directly into 
a Miniati (2001) shows that the kaon channel provides only a small correction to the 
dominant pion channel except at very high secondary energies (where the total number of 
particles is very small regardless). Second, we assume that the ICM is composed entirely 
of hydrogen, while in reality ~ 25% of the mass is contained in helium nuclei. Mannheim 
& Schlickeiser (1994) show that modifying the ambient medium to a normal interstellar 
medium composition increases the tt production rate by < 30%. 

The production rate depends critically on the input flux of cosmic ray protons 

Ulp) = ^/3(7p)iVp(7p), (16) 

where Np^'jp) is the differential number density of cosmic ray protons with Lorentz factors 
between (7^,7^ + rf7p). The protons are expected to be accelerated either by the accretion 
shock surrounding the cluster or by shocks generated during merger events. Such shocks are 
collisionless and non-relativistic, and hence could efficiently accelerate particles to relativistic 
energies (Bell 1978; Blandford & Ostriker 1978; Blandford & Eichler 1987), as observed 
locally in supernova remnants (Koyama et al. 1995, 1997; Tanimori et al. 1998; Muraishi et 
al. 2000; Butt et al. 2001). 

We normalize the proton number density by requiring that the kinetic energy density 
in relativistic protons be a fraction ^cr of the total thermal energy density Wth = Sne/csT of 
the cluster: 

"^pC^ j ilp - ^)Np{lp)d-'ip = CcKUth- (17) 

Shock acceleration models predict that the distribution function of relativistic particles is 
oc p^*, where p is the particle momentum, s = (r + 2)/(r — 1) and r is the shock compression 
ratio (Blandford & Eichler 1987). The differential number density of accelerated protons 
may therefore be written as 

^pilp) — ^P'O _ ]^ys+l)/2 7p,min < 7p < 7p,max, (18) 

where 7p,min and 7p,max are the minimum and maximum Lorentz factor to which a proton can 
be accelerated in the cluster shocks. For strong shocks in a gas with adiabatic index F = 5/3, 
s = 2. Miniati et al. (2001b) argue that the accelerating shocks in cluster environments are 
relatively weak, with Mach numbers < 4-5, corresponding to power law indices s > 2.25. We 
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consider the two cases s = 2 and s = 3 in the following; these should bracket the (uncertain) 
distribution in cluster environments. The normalization constant Np^Q is 



_ 3.2 X 10-^ f^cR 

^^p,0 — 7^ JT'e-'keV 



CCR 



0.1 



cm 



(19) 



where C^^ depends on the power-law index of the cosmic ray distribution function through 
equation (17); for example, 
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+ 



7p,r 



+ 1 



Tp, min + 1 7p,max + 1 



(20) 



(21) 



.(7p,max + l)(7p , mm 1 

Our fiducial value of ^cr = 0.1 is probably a lower limit to the energy deposited in cosmic 
ray protons (Fields et al. 2001; Butt et al. 2001). 

The maximum energy to which protons may be accelerated is relatively well determined. 
Two factors limit 7p,max: the acceleration time, tacc ~ i^lc/v^, must be smaller than both 
the age of the cluster and the cooling timescale; here Vs is the shock velocity (assumed to 
be close to the sound speed in the cluster) and is the proton Larmor radius. The cooling 
time of the protons (which is determined primarily by pion production) exceeds the age of 
the universe, so the first condition is the critical one, yielding 7p, max ~ 
where S^g is the shock magnetic field in microgauss. Note that because 7p,max ^ 1, C3 is 
essentially independent of its precise value, while it enters C2 only logarithmically. 

The value of 7p,min is less certain. Most analytic models of diffusive shock accelera- 
tion assume that the particles retain a Maxwell- Boltzmann distribution below the minimum 
injection momentum pmin and a (broken) power law above that value (e.g., Kang & Jones 
1995; EUison, Berezhko, & Baring 2000). Normally pmin is a few times the mean thermal 
momentum behind the accelerating shock, i.e. ~ 2ci(mpA;BTsh)^^^, where Tgh is the 

postshock temperature. In principle, Ci ~ 2-3 is determined by the assumed energy density 
of the cosmic ray component. However, because pmin is located in the exponential tail of the 
Maxwellian distribution, it depends only very weakly on the cosmic ray energy density. We 
therefore fix Ci ~ 2.3, a value inferred from numerical simulations of the shock acceleration 
process (Gieseler, Jones, & Kang 2000). Then 



7p,i 



3.4 X IQ-'n^v (^) ( — 



■keV 
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where M. is the mean Mach number of the shocks. In our calculations, we fix At = 4, as 
suggested by Miniati et al. (2001b). Note that, particularly in the s — 3 case, a non-negligible 
fraction of the cosmic ray energy may be carried by protons with 7^ < 7th. 

The positron production rate via proton collisions Qppij) is shown in Figure 1 for input 
proton spectra with s — 2 (solid curves) and with s — 3 (dotted curves). Because Qpp is 
strictly proportional to ^cR^e' normalize our results by this quantity. The dependence on 
Tg is slightly more subtle, so we show results for ksT^ = 1, 5, and 10 keV, from bottom to 
top. The ambient temperature does not affect the shape of the positron spectrum (because 
in all cases 7p,min < 7th), but it does affect the overall normalization factor Np^Q [cq. (19)] 
and the fraction of cosmic rays with 7^ > 7th [eq. (22)]. Nevertheless, our results differ only 
slightly from the Qpp oc Tg scaling that one would naively expect from equation (19). Figure 
1 shows that the positrons are produced with a characteristic 7+ ~ 200; this indicates that 
> 99% of the energy deposited in positrons will be lost to coohng before annihilation. At 
high positron energies, Qpp oc 7"*. This occurs because (7^+ is only a weak function of proton 
energy for 7^ > 1. 

For reference, the total positron production rate via proton collisions is Qpp ~ 4.5 x 
10~^^neA^p,o cm~^ for an s = 2 cosmic ray spectrum and Qpp 1.5xlO~-^^neA^p,o cm~^ s~^ 
for an s = 3 spectrum. 

2.3.2. Other Positron Production Mechanisms 

Cosmic ray protons can create positrons through two additional channels: photopair 
production (p + 7 — > p + e"^ + e") and photopion production (p + 7 — > n + tt"^) followed 
by TT"*" decay. In the diffuse IGM, these processes are unimportant except for the highest 
energy cosmic rays because photopair (or photopion) production requires 7p ^ 7 x 10^ (or 
7p ^ 10^^) if the photon originates from the CMB. However, the high-temperature radiation 
field in clusters lowers these thresholds dramatically to 7^ > SllTj^y for photopair production 
and 7p ^ 8 X lO^Tj^y for photopion production. These channels could therefore in principle 
play an important role in positron production in galaxy clusters. The effects of the lowered 
thresholds have not to our knowledge been addressed explicitly in the literature, so we will 
estimate the total production rates here. We assume that the ambient photon field arises 
from bremsstrahlung radiation, with a differential energy density Ui, — Airrcsl^ /3c, where Tc 
is the core radius of the cluster and el^ is the bremsstrahlung spectral emissivity function 
(Rybicki & Lightman 1979). 

For either of these processes, the particle creation rate for a proton traveling with Lorentz 
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factor 7p may be written as (Mannheim & Schlickeiser 1994) 



dx dn /"^^^^ 



27p A' /27p (ia; 



dx' x' (t{x'), 



(23) 



where dn/dx is the (isotropic) photon distribution function, a is the relevant coUision cross- 
section, and x[^ is the minimum photon energy in the proton rest frame needed to create 
the specified particle(s). If we approximate a as a constant and the exponential cutoff in 
the brcmsstrahlung photon spectrum as a step function with cutoff ksTe, then the integrals 
may be easily evaluated. We find that the positron creation rate due to a single proton is 



~~dt 



Kca 



In 



7p 
7* 



+ 



- 1 



(24) 



where X is a constant evaluated below and 7* = x[-^/2kBTe is the minimum proton Lorentz 
factor required to produce the particle(s) given our assumed photon spectrum. 

If we assume a cosmic ray spectrum Np{'jp) = Np^oj^"^, maximizing the effects of protons 
with 7p > 7*, the total rate of positron production Q is 



Q 



KcaNp^o 
37* 



(25) 



For the photoproduction processes we consider, a^e+e- ~ 6 mbarn (Chodorowski, Zdziarski, 
& Sikora 1992), while (7^^+ ~ 150 /xbarn (Bhattacharjee & Sigl 2000). Evaluating K based 
on the assumed energy density of the cluster, we then find 

~ ' X 1°-"^-^- {j^y cm-3 (26) 

and 

Q,,. ~ 5 X 10-^«N,,T^i {^^y (^) cm-3. (27) 

Thus, Q^T^ ^ ^76+6- Qppj and the two photo-collision processes may indeed be 
neglected in even the most luminous clusters. (We find this conclusion to hold regardless of 
the detailed photon spectrum assumed, provided that the total energy density is fixed.) Note 
that, despite our conclusion that Q^^, <^ Qje+e-, photopion production is a more significant 
energy loss mechanism for high-energy protons than is photopair production. The difference 
in significance results from cosmic ray protons with 511 < 7p7kev ^ 8 x 10^ being energetic 
enough to produce pairs (but not pions) and from the greater inelasticity of photopion 
production as compared to photopair production (Mannheim & Schhckeiser 1994). 
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3. Calculation of the Positron Distribution Function 

3.1. Positron Injection by an AGN 

For positrons which are injected impulsively by a short-lived AGN embedded in the 
cluster, the source term Q{'y) = in equation (1). If we further neglect annihilation of 
the cooling positrons, the calculation of the positron distribution function is straightforward 
(Sarazin 1999). This approximation is vahd because for any 7 > 7eq, the annihilation time 
is much longer than the cooling time. (See §4.1 for further discussion of this point.) 

The calculation proceeds as follows. We first assume that the AGN injects positrons into 
the cluster with a distribution function N+{'y,ti), where U is the injection time. Consider a 
positron with Lorentz factor 7j at time U. At a later time t, the positron will have cooled to 
a Lorentz factor 70. Number conservation demands that 

/•OO /"OO 

/ 7V+(7',t)d7'= / N+{^\U)d^'. (28) 

Differentiating, we find (Sarazin 1999) 

7V+(7o,t) = 7V+(7,,t,)^ , (29) 

"70 t 

with 7j computed using the cooling function 6(7, t) given in §2.1. 

At any time t > ti, the number of positrons that have not yet cooled Uuit) may be 
calculated by integrating the distribution function over all Lorentz factors 7 > 7eq. The 
number density of thermalized positrons n+{t) then evolves according to 

where the first term is the rate at which positrons thermalize and the second term is the 
rate at which they annihilate. Thus 



e^(7eq)t' 



dt 



(31) 



We calculate the initial distribution function iV+(7,tj) by expressing the total positron 
energy density (including the rest mass) immediately after injection as 

/L T 
7iV+(7, ti)d'^ = Cagn-^, (32) 
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where Lk is the total mechanical luminosity of the AGN, r is the lifetime of the AGN 
(in a single duty cycle), V is the volume over which the relativistic particles are mixed, 
and ^AGN is the fraction of the injected energy in the positron component (see below). For 
typical parameter values, the energy input from a central AGN is significant compared to the 
thermal energy stored in the cluster core and could, in some cases, balance the energy lost 
to cooling flow radiation (Bohringer et al. 2001b). If we assume that the positrons initially 
have a power law distribution in momentum, 
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the normalization constant n, is 
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where rmix ~ is the radius out to which positrons are distributed and (7^gn (depends on 
the power-law index of the distribution and is determined by equation (32): 
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In equation (33), 7+, max ~ 4 x 10'''Al(-Bs/j^^G^kev)''^^^ is fixed by equating the accelera- 
tion time to the IC coohng time. Bgh thus refers to the magnetic field in the shocks inside 
the accelerating region of the jet system; for simplicity, we take Bgh — B. We note that 
the magnetic field in the accelerating region may have considerably larger values than the 
ambient cluster field because of the powerful magnetic fields in AGN jets (e.g., Begelman, 
Blandford, & Rees 1984). However, the dependence on 7+, max is only logarithmic and as- 
suming a larger B value introduces only a small change to our results. We also assume that 
~ 10, a value appropriate for the hotspot velocity at near the end of the active phase 
[see, for example, the model in Furlanetto & Loeb (2001)]. 

As above, we determine 7+^ min assuming that the minimum cosmic ray momentum is a 
few times the mean thermal momentum. For positrons, this implies that 



<„,-l=0.74r.v(||)^(^) 



(37) 
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Note that (7+,min — 1) is of order unity for reasonable cluster parameters. 

The distribution function in equation (33) should be thought of as referring to the 
positron component once it eventually escapes the entire region of interaction between the 
relativistic jets and the ambient medium. It is not necessarily the same distribution the 
positrons have while inside the jet, where cooling may be efficient (Begelman, Blandford, & 

Rees 1984) but where the positrons are expected to have a bulk Lorentz factor 7jet ~ 10. 
As they travel along the jet and particularly when they pass through the hotspot, surviving 
positrons will presumably be re-accelerated to a fresh distribution given by equation (33). 
The normalization factor ^agn thus represents the net fraction of the interaction energy 
Lkt transferred to the escaping positrons. Note that we assume that all of the positrons 
are accelerated by the shocks. If only a fraction are actually accelerated, the total number 
of positrons (and hence the annihilation signal) would only increase. 

Figure 2 shows iV_|_(7,t) in the AGN scenario at several different time intervals after 
injection. We present results for an s = 2 injection spectrum (sohd curves) and for an 
s = 3 injection spectrum (dotted curves); in each case, curves show to — ti — (1, 2, 3, 5) x 
10^ yr, from top to bottom. Here to is the (present-day) age of the universe in our assumed 
cosmology. All cases assume what we will refer to as our "standard cluster parameters:" 
Lk = 10^5 erg s-\t = 10^ yr, = IQ-^ cm-'\ keT^ = I keV, B ^ 3 fiG, r^ix = 200 kpc, 
and ^AGN = 0.1. Note that N^(^,t) oc ^agn Lkti"^^^^, with a more subtle dependence on 
He (see below). This scaling holds as long as our approximation of instantaneous positron 
injection is valid, namely for r < 3 x 10® yr. The dependence on temperature enters only 
through 7+,mm- 

The main feature in the distribution function is the sharp cutoff at a Lorentz factor 
7max(^)- At any given time after injection, this corresponds to 

7max(^) = lim 7(7j,t,tj), (38) 

where '-f{'ji,t) is the Lorentz factor at time t of a positron injected with 7^ at time ti. The 
limit is well-defined because 6(7) oc 7^ for large Lorentz factors; however, we do not solve 
for it in closed form here. Sarazin (1999) gives expressions for 7max(^) in some simple limits. 

Figure 3 shows the effects of varying the cluster core density rie on the distribution 
function of the uncoolcd positrons. Again, the solid (dotted) curves show results for an 
s = 2 (3) injection spectrum. All assume to — t^ = 10^ yr and, aside from the density, 
standard cluster parameters. From top to bottom, the curves in each case assume Ug — 
10~^, 10~^, and 10~^ cm~^. Perhaps surprisingly, N^{'y) decreases as Ue increases. This 
non-intuitive behavior is a result of the cooling rates. Above the coohng "bottleneck" at 
7 ~ 200(ne/10~^ cm~^)^/^, the cooling time is independent of density [see equation (5)]; 
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below this Lorentz factor, the coohng time is inversely proportional to density [see equation 
(7)]. In other words, once positrons enter the Coulomb cooling regime, they thermalize 
at a rate proportional to fig. Thus, near the bottleneck iV+(7) is independent of density, 
while below it A^+(7) actually decreases with increasing density because positrons enter the 
thermalized component faster. (This density dependence does not hold for the thermalized 
positron density n+.) 

The evolution of function of cosmic time is shown in the bottom panels of Figures 

4 and 5. Figure 4 shows results for an s = 2 positron injection spectrum; the solid curves 
assume Ue = 10~^ cm~^ with t^ — ti = 1, 2, 3, and 5 x 10^ yr, while the dotted curve assumes 
He = 10~^ cm~^ with tQ — ti = 5 X 10^ yr and the dashed curve assumes Ue = 10~^ cm~^ with 
to — ti = 5 X 10^ yr. Figure 5 shows a similar set of results for an s = 3 injection spectrum. 
In all cases, we assume standard cluster parameters (except for the electron density). We 
find that steepening the power law slope s increases the number of thermalized positrons 
because the total number of positrons produced by the AGN increases as s increases. 

At very early times after injection, the thermalization rate exceeds the annihilation rate 
and n+ rises rapidly, particularly for a steep injection spectrum in which the number density 
of particles is heavily weighted toward small 7. However, once the low-7 positrons cool, the 

thermalization rate falls rapidly with time (because fewer positrons begin with large Lorentz 
factors) and annihilation begins to dominate. Therefore, at late times n+ oc e~^^'^'"^^^^~^'-\ 

The dependence on density is more complex. The total number density of injected 
positrons (and hence the maximum value of n+) is independent of Ug, and so is the peak value 

of n+. However, because fecoui oc Ug, the value of n+ peaks later for smaller ambient densities. 
Furthermore, the annihilation timescale is Tann ~ ^(7cq)^^ ~ 4 x 10'^(ne/10^"^ cm^"^)^-*^ yr; 
for rig = 10""^ cm~^, Tann ^ and n+ does not change substantially over the lifetime of 
the cluster. 

Note that the distribution functions derived in this section depend weakly on the cos- 
mological redshift through the corresponding dependence of the IC cooling rate. 

3.2. Steady Injection by Cosmic Rays 

In the case of steady injection, with Qppij) computed as in §2.3.1, we can calculate the 
positron distribution function reached at steady-state (i.e., when the total prodiiction rate 
of positrons equals the annihilation rate). In this case, all quantities in equation (1) are 
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independent of time and the distribution function of uncooled positrons is 



1 r°° 



QiY) 



(39) 



6(7') ' 



where 



^(7) = exp 




di \dh 



A{i) 



(40) 



h{i) [di 



and 7i is an arbitrary constant. Note that, if we neglect annihilation, this reduces to the equi- 
librium distribution found by Sarazin (1999). The number density of thermalized positrons 
may then be found by requiring the total number density of positrons to be independent of 



Our scheme requires that enough time has elapsed since the cosmic ray proton compo- 
nent first formed for the annihilations to reach equilibrium with production. Because the 
annihilation cross section decreases rapidly with 7_|_, this essentially requires that the ther- 
malized component has grown enough that A(7cq)n+ ~ Qpp. This requires a time Tcq ~ Tann- 
We therefore expect equilibrium to be a satisfactory assumption for ng > 3 x 10^^ cm~^, a 
condition well-statified by rich clusters. 

Figure 6 shows the equilibrium positron distribution A^+(7) with QppIt) from §2.3.1. 
We show results for an s = 2 cosmic ray proton spectrum (solid curves) and for an s = 3 
proton spectrum (dotted curves). In each case, curves correspond to /c^Tg = 1, 5, and 
10 keV, from bottom to top. All assume Ue — 10~^ cm~^, B — ^ fxG, and .^cr — 0.1. The 
distribution function is strictly proportional to ^cr- Although Tg does affect 7p,min and hence 
the fraction of cosmic ray protons able to produce pions, we find that N^{'y) is still very 
nearly proportional to Tg. 

Figure 7 illustrates how A^+(7) depends on the density of the ambient medium. As 
before, the solid (dotted) curves show results for an s = 2 (3) proton spectrum. The cluster 
is assumed to have ^cr — 0.1, .6 = 3 //G, and ksTg — 10 keV, with — 10~^, 10~^, and 
10~^ cm~^, from bottom to top. We see that the peak in the distribution function moves 
to larger Lorentz factors as density increases and Coulomb cooling becomes more efficient. 
Note, however, that the cluster with rig — 10~^ cm~^ may not yet have reached steady-state. 

In all cases, A^+(7) cx: 7 at high energies where synchrotron/IC cooling dominates 
and annihilation can be neglected (Sarazin 1999). The distribution function peaks where 
the efficiency of both synchrotron/lC and Coulomb cooling is low (7 ~ 100) and then falls 
again at small Lorentz factors as Coulomb cooling takes over. 



time: 




7eq 



(41) 
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The equilibrium density of thermalized positrons n+ is not shown in Figures 6 and 7 
because it scales simply with the cluster parameters: 

^ 8 X 10-"X«r,.v (35^) (If) cn.-3, (42, 
where = 1 (2.5) for an s = 2 (s = 3) proton spectrum. 

4. The Annihilation Spectrum 

In the following subsections we calculate the emissivity of annihilation photons from the 
ICM using equations (10) and (13) for both the AGN (§4.1) and the steady injection (§4.2) 
scenarios. In each case, we use the distribution function iV_|_(7) found in the appropriate 
subsection of §3. 



4.1. Positron Injection by an AGN 

For injection by an AGN, we ignored annihilation of non-thcrmalized positrons when 
calculating iV+(7,t), although we did include annihilation of the thermalized component. 
Here we calculate the instantaneous emissivity expected for the resulting distribution. 

Figure 8 shows the annihilation emissivity {dn^ / de-ydt) for the same cases as in Figure 
2. The continuum, in which the emissivity is approximately oc e~^, is generated by the 
uncooled positrons, while the line is generated by the thermalized positrons. The sharp 
cutoffs in the continuum occur because of our cosmic ray approximation in which 7_ = 1. In 
this approximation, only photons with [7+,max(l - /3+,max) + l]/2 < e^jm^c^ < [7+,max(l + 
/3+,max) + l]/2 can be generated (Svensson 1982); in reality, the non-zero high-energy tails 
of the ambient electron distribution (and collisions with relativistic cosmic ray electrons) 
will cause the cutoffs to be somewhat smoother. The time evolution of line emissivity is 
discussed in detail below. 

Figure 8 justifies our approach of neglecting annihilation in calculating the distribution 
function. Consider, for example, the case in which to — U — 10^ yr for an s = 2 injection 
spectrum. The total annihilation emissivity in the continuum is ~ 5 x 10~^^ cm~^ s~^, 
yielding a total number of annihilations of non-thcrmalized positrons since injection of ~ 
8 X 10~^^ cm^^. We calculate the total density of non-thermalized positrons to be ~ 2 x 
10~^° cm~^, indicating that only ~ 5% of the positrons annihilate before cooling. 

The line emissivity at a given time is proportional to ^AGNLxTr^y., and so the total 
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luminosity of the cluster is independent of rmix- The dependence on ambient density is 
somewhat more comphcated. This is displayed in Figure 9, which shows the annihilation 
emissivity for clusters with densities Ue = 10"'^ cm~^ and rig = 10~^ cm~^ for s = 2 and 
s = 3 positron input spectra (solid and dotted curves, respectively). Standard cluster pa- 
rameters (with the exception of density) are assumed, with to — ti = 10^ yr. The continuum 
emissivity depends on the shape of N^{j,t), which has a non- negligible dependence on Ug. 
Note that more continuum photons are produced in lower density media. This is because, as 
shown in Figure 3, a larger fraction of the positrons have thermalized by this epoch in high 
density clusters, and the thermalized component contributes only to the line emission. The 
dependence of the peak line luminosity on density is very different from that of the contin- 
uum (see below). The cluster temperature affects both the normalization of the continuum 
(although only relatively weakly) and the width of the annihilation line [the full width at half- 
maximum is approximately ~ 32T^/,y keV; see equation (13)]. This temperature-dependent 
line broadening is displayed graphically in Figure 10 below. 

Observationally, the critical observable is the emissivity of the annihilation line; this 
quantity varies strongly as a function of the source age. We show the evolution of the line 
emissivity nunc with time in the top panels of Figures 4 and 5 for the same cases shown in 
the bottom panels. The emissivity is proportional to the number of thermalized positrons, 
so at late times hune oc nee"^^'''''^)^*"**). Thus the peak signal decreases as density decreases, 
but it also fades with a timescale Tann oc 

The cmissivities calculated in this section depend weakly on redshift through the IC 
cooling rate. 



4.2. Steady Injection by Cosmic Rays 

In the steady injection scenario, we have calculated A^+(7) self-consistcntly, including 
annihilation. The resultant photon cmissivities arc shown in Figures 10 and 11 for s = 2 
and s = 3 proton spectra, respectively. In each case, the solid curves show the annihilation 
emissivity, with ksTe = 1 keV (bottom) and ksTe = 10 keV (top). 

In these figures, we also show the inverse Compton emissivity determined from the 
equilibrium positron spectrum A^_|_(7). In the energy range of interest (7 ~ 10^), A^+(7) has 
assumed its asymptotic power law form, and the inverse Compton spectrum can be estimated 
easily (Rybicki & Lightman 1979): 

- 5 X 10 /.T,ev J (iF^^) ) ' ^ ^ 



de^dt 
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where Is — I (10^) for an s = 2 (s = 3) proton spectrum. This is shown in Figures 10 and 11 
by the dotted hnes. Clearly, the IC emission from the positrons themselves is a substantial 
contaminant, and in fact hides the annihilation radiation completely for relatively shallow 
proton spectra. Thus, only if the protons are accelerated by very weak shocks (so that s ~ 3) 
is the annihilation radiation from the secondary positrons observable even in principle. 

In the steady injection case, [dn^/de^dt) oc Np^QUe oc ^cR^e- ^^^^ ambient density also 
determines the energy at which the continuum begins to decline rapidly; the break is located 
at ^p]^mf.c^, where is the location of the peak of A^+(7) and is determined by the 

efficiency of Coulomb cooling. The temperature affects the normalization of the continuum 
and line (nune oc Tg) as well as the determination of the width of the Une component, as 
described in the previous subsection. 



5. Discussion 

We have calculated the signals expected from positron annihilation in galaxy clusters 
for positrons injected as primaries by embedded radio jets and those produced as secondaries 
in collisions between cosmic rays and thermal protons in the cluster. The former case is of 
interest in constraining the matter content of relativistic jets, while the latter is of interest in 
measuring the cosmic ray content of the ICM. The positron annihilation line is particularly 
interesting for these purposes because typical cluster temperatures /c^Te ~ 1-10 keV are 
just in the range in which positrons annihilate efficiently without forming positronium. We 
calculate the annihilation rates assuming that the cooling positrons mix efficiently with the 
ambient medium on cosmological timescales. 

For positron injection by a single AGN, the peak emissivity in the annihilation line is 

nHne,sing.e ~ 2X10 X, [^^^^^3 ) ) [^^^ ) 

X ( in4.^'^ -i ) ( T7^ ) s-^ (single AGN event) (44) 



10^5 erg s-V VlO^ yr 

with X^'^^ = 1 (10) for an s = 2 (3) positron injection spectrum. The dependence on s 
appears primarily because, with the normalization procedure described in §3.1, the total 
number of e'^ produced increases as s increases. We note that rich clusters can have core 
densities rig ~ 10~^ cm~^ and ^agn may reach a value of several tens of percent for jets with 
a low bulk Lorentz factor. The time lag between injection and peak annihilation depends 
principally on the electron density rie in the cluster core: the characteristic annihilation time 
of the positron population is Tann ~ 4 x 10^(ne/10~^ cm~^)~^ yr, and the emissivity maintains 
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the peak level for roughly this time period (see Figures 4 and 5). Therefore, although dense 
clusters produce the strongest signals, they fade relatively quickly. 

The longevity of the positron population in typical clusters indicates that the signal 
can be enhanced if we consider either multiple injection epochs or multiple AGN in a single 
cluster. A model of the former type has been recently advocated by Bohringer et al. (2001b) 
as a way of balancing the cooling flow radiation in some clusters. These authors suggest 
quasi-periodic mechanical energy injection from the central galaxy throughout the lifetime 
of the cluster, with injection epochs lasting ~ 10*^ yr occurring every Atinj ~ 10^ yr. In the 
limit in which Tann ^ Atinj, the annihilation cmissivity will reach a quasi-steady state with 
*^iine ~ ''7 ^^^line, single) whcrc Tj = Atinj/Tann IS the pcriod betwccu outbursts in units of the 
annihilation time of the cluster. Here and r are to be interpreted as the mechanical 
luminosity and lifetime of a single outburst. 

An even more interesting possibility is positron injection by galaxies throughout the 
cluster. While the vast majority of these galaxies are not active at the present day, recent 
studies have found evidence for relic supermassive black holes in nearly all bulge-dominated 
galaxies (Magorrian et al. 1998; Gebhardt et al. 2000). This, together with modeling of the 
quasar luminosity function, suggests that nearly all galaxies once hosted a quasar (Haiman 
& Loeb 1998; Hachnclt & Kauffmann 2000). Current observations give the relation Mbh = 
7.8 X 10^(LB,buigc/10^° Lb0)1-°^ Mq (Kormendy & Gebhardt 2001). Because the observed 
black hole mass-bulge luminosity relation is nearly linear, for the purposes of an estimate it 
suffices to assume a linear relation and scale our results with the total cluster core luminosity 
Lb, ci- A typical rich cluster core has a total i?-band luminosity of L^^ci ~ 10^^ L^© (Peebles 
1993). [This luminosity corresponds to ~ lOOL^ galaxies in the cluster, where L*^ is the 
characteristic luminosity of galaxies in the Schechter function (Yasuda et al. 2001).] We 
further assume a mass-to-energy conversion efficiency £bh during the black hole formation 
process and that a fraction fx of this energy is released in outflows. We expect £bh ~ 0.1, 
and observations indicate that Jk ~ 0.1 (Hooper et al. 1995; Furlanetto & Loeb 2001, and 
references therein). Therefore, if we make the extreme assumption that all injection events 
occur simultaneously, the peak line emissivity would be 

nnne ~ 9 X 10 X, [j^z^^) [-^ ) [^^^ ) 

^ / — B,ci — \ g-1 (multiple simultaneous AGN) (45) 

VO.l 0.1 10^2 Lbq J ^ \ f / V ; 

with X^^^ deflned as above. Of course, we must keep in mind that the quasar era peaked 
at z ~ 2 (Pei 1995; in our cosmology, to — ti ~ 10^° yr) so a substantial fraction of the 
population may have annihilated before the present day. In the opposite limit, in which the 
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source evolution time is much larger than the annihilation time, but in which the annihilation 
time is in turn much larger than the time between injection events (or Hq^ ^ Tann 3> Atinj), 
a quasi-steady state will be reached as described in the previous paragraph, with a steady 
annihilation emissivity of approximately r]~^ times that of a typical AGN in the cluster. 

For positron production by cosmic rays, the steady-state emissivity in the annihilation 
line is 

'T'line ~ 10~^°^?^ 2kev f — cm~^ s"-*-, (cosmic rav Secondaries) (46) 

\ 10"'^ cm V 0.1 / 

with X^^ = 1 (2.5) for an s = 2 (s = 3) proton spectrum. For an s = 2 input spectrum, the 
energy lost to positron annihilation over the age of the universe is ~ 10~^C/cR(^e/10~^ cm~^), 
where Ucr is the cosmic ray energy of the cluster. The emission mechanism is inefficient 
because < 10~^ of the cosmic ray energy goes into the positrons and because the characteristic 
initial energy of the positrons is » mgC^ (see Figure 1) so that most of the initial positron 
energy is lost to cooling radiation before annihilation occurs. We therefore see that direct 
positron injection by even a weak AGN {Lk ^ 10'^^ erg s~^) would overwhelm the signal 
from secondary positrons produced by cosmic ray protons. 

The observability of the line is of course the critical question. The flux at earth from a 
cluster with a positron mixing radius r^i^ is 

= 8 X 10-^ ( -^^X ( .R. 1 ' i r. s 1 ^ photons cm'^ s-\ (47) 

V200 kpc/ VlOOMpc/ Vl0~^cm-3s-V 

where D is the luminosity distance to the cluster. In the AGN injection case, the emissivity 
is oc r~fx, and so the flux is independent of the mixing scale (rather it depends on the total 
number of positrons injected by the AGN). For steady production by cosmic ray protons, the 
emissivity is independent of the volume of the cluster core and the flux scales in proportion 
to the mixing volume. 

Before discussing the prospects for detection of this signal with upcoming instruments, 
we must consider possible contaminating backgrounds. First, as described in the previous 
subsection, the relativistic positrons (and electrons) with 7 > 10^ in the cluster generate IC 
radiation at 511 keV. However, cooling depletes this population rapidly; after a time interval 
t — ti^ 1.5 X 10^ yr, the maximum Lorentz factor 7max < 10^ and the electrons and positrons 
can no longer produce IC radiation at 511 keV. However, in the case of steady injection from 
cosmic ray protons, the IC background originates from the equilibrium distribution itself. 
Figures 10 and 11 show that in this case the IC background can overwhelm the annihilation 
line unless the protons have a sufficiently steep {s ~ 3) injection spectrum. 
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Second, we must also consider the diffuse extragalactic background, for which the flux 
at ?a 511 kcV is Fbkgd ~ 2 x 10~^ photons cm~^ s"^ keV~^ sr~^ (Watanabe et al. 1997). 
For a cluster of a fixed size, the ratio of the line flux from the cluster to the background flux 
is independent of cluster distance so long as ^ <^ 1, 



Thus, the expected signal from AGN-injected positrons is well above the diffuse background. 
However, the signal from secondary positrons generated by cosmic rays will be hidden by 
the background. 

While the prospects for observing the annihilation of secondary positrons in the foresee- 
able future are small (even if the annihilation line can be observed over the IC background), 
positrons injected by AGN may soon be detectable with spaceborne instruments. The IN- 
TEGRAL satellite^, expected to be launched in October 2002, will have spectral capabilities 
in the energy range of interest. The SPI instrument is expected to have a 3a line sen- 
sitivity ~ 5.1 X 10~® photons cm^^ s^^ given an integration time of lO'^sec, but its poor 
angular resolution (2.5°) may lead to background contamination. The IBIS instrument, 
with 12' resolution, is better suited to cluster detection, but it has a 3a line sensitivity of 
only ~ 2 X 10"^ photons cm~^ s~^ (again for an integration time of 10^ sec). An even more 
powerful search could be conducted with EXIST^, a proposed all-sky hard X-ray survey mis- 
sion. It has an expected 5a line sensitivity of ~ 5 x 10~^ photons cm~^ s~^ in the relevant 
energy range (assuming an integration time of 10^ sec, the mean exposure time planned for 
any point on the sky in the mission), and it has an excellent angular resolution of 5'. If 
the positrons are injected by AGN, these sensitivity limits are close to the signal we predict 
for nearby (D < 100 Mpc) clusters with powerful AGN, multiple injection epochs/galaxies, 
or steep injection spectra. Deep exposures with INTEGRAL or statistical analyses taking 
advantage of the full sky coverage of EXIST may reveal weaker sources as well. 

A particularly interesting source is the nearby Virgo cluster, at a distance ~ 20 Mpc. In 
the AGN scenario, the annihilation signal in Virgo is comparable to the detection threshold 
of all three instruments listed above. The expected peak flux is 




(48) 





photons cm ^ s 



-1 



(49) 



^See http://astro.estec.esa.nl/SA-general/Projects/Integral/integral.htnil 
^See http://exist.gsfc.nasa.gov 
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where we have used fiducial values for the luminosity of M87 estimated by Bohringer et al. 
(2001b) and for rif. by Nulsen & Bohringer (1995). Here 'q~^ represents the contribution from 
past AGN phases of cluster galaxies [see the discussion accompanying equation (45)]; in the 
best case, it could represent an enhancement of more than an order of magnitude. If rmix 
is large, the high resolution instruments may even be able to map spatial variations in the 
positron component (and indicate whether the positrons are injected solely by the central 
galaxy or by a larger number of galaxies during the quasar era). 

Continuum 7-ray emission from M87, the dominant galaxy in Virgo, may contaminate 
the annihilation line signal. To date there have been no observations of Virgo in the relevant 
energy range. To estimate the contamination from M87, we assume that observed power- 
law spectrum in the 2-10 keV range extends to 511 keV; using the recent observations of 
Bohringer et al. (2001a), we expect M87 to produce a fiux in the spectral regime of the 
annihilation fine of ~ 4 x 10~^T]^/y photons cm~^ s~^ (the temperature of the core of Virgo 
varies with radius between /csT'e = 1-3 keV; see Bohringer et al. 2001a). Thus, even if 
the continuum emission from M87 cannot be removed through the use of a high-angular 
resolution instrument, the annihilation line should still be visible. Of course, because M87 is 
still active, any positrons produced in the current outburst phase have most likely not had 
sufficient time to cool. We would therefore expect annihilation line emission only if there is a 
relic population of positrons either from earfier outbursts of M87 or from the other galaxies 
in Virgo. 

Another interesting source is Centaurus A, the nearest {D k, 3.5 Mpc) bright radio 
galaxy to the Milky Way. The radio structure is one of the largest in both apparent and 
absolute size, covering an area 8° x 4° on the sky (or 480 kpc x 240 kpc) and is therefore a 
very attractive candidate for observations (Israel 1998). The expected fiux in the positron 
annihilation line is 



where we have scaled Lk to the approximate bolomctric luminosity of the Centaurus A nu- 
clear source (Chiaberge, Capetti, & Celotti 2001) and r to the minimum source age estimated 
by Saxton, Sutherland, & Bicknell (2001). Recent data from the Chandra X-ray Observa- 
tory indicates that the gas density around Centaurus A is rig ~ 10^^(r/5 kpc)^^"^'^ cm~^ 
at 0.5 kpc <^ r < 10 kpc (Kraft et al. 2002). The signal is potentially stronger than that 
of M87; however, several caveats are in order. First, the gas near the central source has 
a temperature ksTf. ~ 0.275 ± 0.03 keV (Kraft et al. 2002), within the regime in which 
positronium formation begins to play a significant role. More importantly, the scale of the 
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raxlio emission suggests that a large fraction of the positrons escape the interstellar medium 
of Centaurus A and mix with the intragroup medium at r > 100 kpc. The inferred virial 
temperature of the group is only ~ 0.07 keV (based on the observed velocity dispersion of 
group galaxies; van den Bergh 2000). At this temperature, the positronium formation rate is 
approximately equal to the free annihilation rate (Crannell, Joyce, Ramaty, & Werntz 1976). 
The gas density will also be much lower at large radii (if the observed power-low decline at 
r ~ 10 kpc continues, then rig ~ 2 x 10~^ cm~^ at r ~ 100 kpc). Therefore it is unclear 
whether annihilation in this region can produce a strong line. In addition, because the source 
is still active, the positrons may not yet have had time to cool or mix sufficiently with the 
surrounding gas. Finally, because the Centaurus A group contains only ~ 30 galaxies (van 
den Bergh 2000), additional positron enrichment from other AGNs is likely to be small. 

In summary, we have shown that, although there are a variety of mechanisms for pro- 
ducing positrons in galaxy clusters, only direct injection by AGNs is likely to produce an 
observable signal. Therefore, we argue that a positive detection of positron annihilation lines 
from clusters suspected of harboring dormant AGNs would be a robust indication that radio 
jets contain an e^e~ pair plasma. 
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Fig. 1. — Normalized positron production spectrum from proton-proton collisions, 
Qpp(7)/(^eCcR)- Shown are source functions for an s = 2 {solid curves) and an s = 3 
cosmic ray proton spectrum {dotted curves). In each case the curves assume ksTg = 1, 5, 
and 10 keV, from bottom to top. 
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Fig. 2. — The distribution function of non-thermalized positrons N+{'y,t) for injection via 
an AGN at various times. Shown are an s = 2 injection spectrum {solid curves) and an 
s = 3 injection spectrum (dotted curves). In each case, curves correspond to At = to — ti = 
(1, 2, 3, 5) X 10^ yr, from top to bottom. All cases assume standard cluster parameters (see 
text) . 
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Fig. 3. — The distribution function of non-thermalized positrons N+{'y,t) for injection via 
an AGN at various ambient densities. Shown are an s = 2 injection spectrum (solid curves) 
and an s = 3 injection spectrum [dotted curves). In each case, curves correspond to rie = 
10"*^, 10"'^, and 10^^ cm''^, from bottom to top. All cases assume to — ti = 10^ yr and 
standard cluster parameters (with the exception of rig; see text). 
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Fig. 4. — The evolution of the thermahzed positron density n+ (bottom panel) and of the 
emissivity in the annihilation line riune {top panel) as a function of cosmic time for an s = 2 
positron injection spectrum in the AGN case. Shown are results for rig = 10~^ cm~^ with 
to — ti — 1, 2, 3, and 5 x 10^ yr from left to right {solid curves), — 10~^ cm~^ with 
^0 ~ = 5 X 10^ yr {dotted curve), and rig = 10~^ cm~^ with — ti — b x IQ^ jr {dashed 
curve). All curves assume standard cluster parameters (with the exception of rig; see text). 
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Fig. 5. — Same as Figure 4, but with an s = 3 positron injection spectrum. 
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Fig. 6. — The equilibrium distribution function of non-thermalized positrons N+i'j) assuming 
steady Qppij), for various cluster temperatures. Shown are an s = 2 proton spectrum {solid 
curves) and an s = 3 proton spectrum {dotted curves). In each case, curves correspond to 
ksTe = 1,5, and 10 keV, from bottom to top. All curves assume — 10~^ cm~^, B — 3 /iG, 
and {cR = 0.1. 
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Fig. 7. — The equilibrium distribution function of non-thermalized positrons N+i'j) assuming 
steady Qpp{j), for various cluster densities. Shown are an s = 2 proton spectrum {solid 
curves) and an s = 3 proton spectrum {dotted curves). In each case, curves correspond to 
He = 10^'^, 10^'^, and 10^^ cm~^, from bottom to top. All curves assume ksTg — 10 keV, 
B = 3 fiG, and ^cr = 0.1. 
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Fig. 8. — The annihilation emissivity {dn^/de^ dt) for injection via an AGN at various times. 
Shown are an s = 2 injection spectrum {solid curves) and an s = 3 injection spectrum 
{dotted curves). In each case, curves correspond to /St — — ti — (1, 2, 3, 5) x 10^ yr, from 
top to bottom in continuum leveL All curves assume standard cluster parameters (see text). 
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Fig. 9. — The annihilation emissivity {dn-y/ de-y dt) for injection via an AGN at various am- 
bient densities. Shown are an s = 2 injection spectrum {solid curves) and an s = 3 injection 
spectrum [dotted curves). In each case, we show results for rig = 10~^ and 10^^ cm^''. 
All curves assume standard cluster parameters (with the exception of rie, see text) and 
to-ti = 10^ yr. 
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Fig. 10. — The spectral emissivity {dn-y / de-y dt) assuming steady Qpp{'^) for an s = 2 proton 
spcctiTim. Solid curves show the annihilation emissivity for a clusters at ksT^ — 1 keV 
(bottom) and ksTe = 10 keV (top). Dotted curves show the inverse Compton emissivity 
determined from A^+(7) for the same temperatures. All curves assume rig = 10~^ cm^^, 
B = 3 nG, and ^cr = 0.1. 
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Fig. 11. — Same as Figure 10, but with an s = 3 proton injection spectrum. 



